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ABSTRACT 

We use the Sloan Digital Sky Survey (SDSS) spectroscopic sample to constrain the projected radial 
distribution of satellites around isolated ~ L* galaxies. We employ mock galaxy catalogs derived from 
high-resolution cosmological simulations to investigate the effects of interloper contamination and show 
that interlopers significantly bias the estimated slope of the projected radial distribution of satellites. 
We also show that the distribution of interlopers around galaxies is expected to be non-uniform in 
velocity space because galaxies are clustered and reside in crowded environments. Successful methods 
of interloper contamination correction should therefore take into account environments of the host 
galaxies. Two such new methods are presented and the most reliable of them is used to correct for 
interloper contamination in analyses of the SDSS galaxy sample. The best fit power-law slope of 
the interloper-corrected surface density distribution of satellites, E(i?) oc R a , in the volume-limited 
SDSS sample is a ~ —1.7 ± 0.1, independent of the galaxy and satellite luminosities. Comparison 
with ACDM simulations shows that the radial distribution of the SDSS satellites is more concentrated 
than that of subhalos around galaxy-sized halos, especially at R < 100/iT 1 kpc. The predicted dark 
matter radial distribution is somewhat more concentrated than the profile of the SDSS satellites, but 
the difference is not statistically significant for our sample. 

Subject headings: cosmology: theory - dark matter - galaxies: formation - galaxies: structure - 
galaxies: fundamental parameters 


1. INTRODUCTION 

In the Cold Dark Matter (CDM) paradigm, satel¬ 
lite galaxies are expected to be associated with the 
dark matter subhalos - halos which lie within the 
virial radius of a larger halo - ubiquitous in the cos¬ 
mological CDM simulations. The abundance and ra¬ 
dial distribution of satellite galaxies can therefore serve 
as a useful test of CDM galaxy formation models, 
constraining the relation between galaxies and subha¬ 
los. In addition, satellite dynamics can provide use¬ 
ful constraints on the total mass distribution in galac¬ 
tic h al os (e.g., IZaritskv fc Whim [1991 IZaritsk v et al l 
jT97LJPrada et_al| 1200.1 Ivan den Bosch et alJ 12004 
Conroy et al 120051 . This, however, requires a good un¬ 
derstanding of how the spatial distribution and kinemat¬ 
ics of satellites and dark matter are related. 

Many recent studies based on numerical simulations 
have shown that the radial distribution of subhalos 
in cluster-sized systems is less concentrated than that 
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of dark matter in the inner sa 20 — 50% of the virial 
radius of host halos, but approxima tely follows the dark 
matter dist r ibutio n at larger rad ii dChigna et al] 1199a 


Colin et all 11999) Ghigna et al 

200(J Suringel et alJ 

20011 De Lucia et al. 20(3 

Diema.nd et al. 2004 

Gao et al. 2004; iNaerai & Kravtsov 20051). Theoretical 


predictions for galaxy distributions in clusters have also 
been accompanied by rapidly improving observational 
measurements (e.g ., HdnetalJ 12004 |Haiia£net_alJl2QQ5l: 
Colhster fc LahavT 120051: lYamTetal.l 120051 : ICoiletal] 
2006), which also find concentrations of galaxy radial 
profiles lower than the concentrations expected for the 
matter distribution of their parent halos. 

The observed distribution of satellite galaxies in 
galactic halos has been studied less extensively. The 
Local Group dwarf population is more radially con¬ 
centrated th an subhalos in dissipationless numerical 
simulations llKravtsov et all 1200414 iTavlor et alJ 12004 
IWillman et al,ll2004) .' a bias that is likely related to the 
phy sics of the formation of the smallest dwarf galax¬ 
ies llKravtsov et alJ 1200414 IPiemand et alJ 120051) . The 
known population of the Local Group satellites is, how¬ 
ever, quite smal l compared to the expected population o f 
CDM subhalos llKlvpin et alJll999bflMoore et al .111 9901) . 
Moreover, the strong radial bias exhibited by the faint 
Milky Way satellites is not expected to apply to the 
brighter satellites (such as, for example, the Magellanic 
Clouds). 

More accurate, statistical constraints on the satel¬ 
lite distribution can be obtained by using galaxy 
redsliift surveys. Several early studies attempted 
to constrain the small-scale galaxy correlation func¬ 
tion by estimating the surface density of objects pro¬ 
jected near galaxies, E(i?) oc R a , finding slopes 
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ranging from a = —0.5 to — 1.25 |lLake_^TxeiIiaind 
l98Gl IPhillipps fc Shank j|1987t IVader fc Sandagdll99H: 

Lorrimer et alJ 11994 ISmith et al.l 12004 iMadore et alJ 
2004 . 

Recently, the availability of large galaxy redshift sur¬ 
veys has allowed construction of large statistical samples 
of parent galaxies and satellites with well defined selec¬ 
tion criteria. The large sample sizes and redshift infor¬ 
mation make it possible to understand the biases and 
completeness of the sample. In addition, isolation crite¬ 
ria for the primaries can be introduced in order to reduce 
the interloper co ntamination and simp l ify th e interpre¬ 
tation of results. Ivan den Bosch et~aL (j2005l l use mock 
galaxy redshift samples derived from large cosmological 
simulations to develop an iterative method of interloper 
rejection for the Two Degree Field Galaxy Redshift Sur¬ 
vey (2dFGRS) and find that the data is generally con¬ 
sistent with the dark matter profile at large projected 
radii, but conclude that incompleteness of close pairs in 
the surve y prevent strong c onstra ints. In an independent 
analysis, ISales fc Lambasl (2005) account for the close- 
pair bias in the data by estimating completeness with 
control samples of objects that are not physically bound 
to the primaries. They estimate the power-law slope of 
the satellite distribution to be a = —0.96 ± 0.03 for pro¬ 
jected radii between 20 and 500 h~ x kpc with a significant 
dependence on morphological type of the parent galaxies 
(a ~ —1.1 for the early type, and ~ —0.7 for the late 
type galaxies). Note, however, that these values of a are 
obtained without any correction for interlopers. 

Given that the satellite distribution can be directly 
probed only in projection, with only limited information 
about positions of likely satellites in three dimensions, 
one has to worry about contamination by interlopers, 
the objects that are not true satellites but are simply 
close to the parent due to projection. Unfortunately, in 
practice it is often tricky to estimate and correct for the 
interloper contamination. This is especially difficult if 
the redshift information is absent as was the case in the 
earliest studies of the satellite distribution. However, 
even in studies in which redshift informatio n is available, 
the in terlopers are often neglected (e.g. . ISales fc Lambasl 
12005 1. Nevertheless, as we show in this paper, the effect 
of interlopers must be corrected for in order to obtain an 
unbiased measurement of the satellite projected radial 
distribution. 

The simplest assumption one can make is that the sur¬ 
face density of interlopers is uniform. The interloper 
contamination can then be estimated by sampling the 
environments around random points in the field. This 
method thus presumes that the volume around a ran¬ 
dom point on the sky and in redshift space contains a 
representative density of interlopers. However, bright 
galaxies are strongly correlated in space and thus can 
be expected to be preferentially located in crowded en¬ 
vironments. One may suspect, then, that the random 
points method can underestimate the interloper number 
density around real galaxies. Therefore, more sophisti¬ 
cated methods, which sample interlopers in the environ¬ 
ments similar to those of the primary galaxies, need to 
be developed. 

In this study, we develop two new methods to estimate 
the contribution of interlopers to the surface density of 
satellites, which take into account the clustering of par¬ 


ent galaxies. We use cosmological simulations to test 
different methods of interloper subtraction and present 
detailed discussion of their strengths and weaknesses. 
We show that interloper contamination can significantly 
bias measurements of the projected radial distribution 
of satellite galaxies. Proper interloper subtraction is, 
therefore, a must in studies of the radial distribution of 
satellites. We use the Sloan Digital Sky Survey (SDSS) 
spectroscopic sample to measure the projected radial dis¬ 
tribution of satellites around nearby bright galaxies, cor¬ 
rected for interlopers. We compare the result to the pre¬ 
dictions for the dark matter and subhalo distribution in 
the ACDM cosmology. 

The paper is organized as follows. In § [21 we discuss 
the interloper contamination and different methods of 
interloper subtraction, testing each of them using mock 
satellite catalogs derived from cosmological simulations. 
We then describe our SDSS spectroscopic galaxy sample 
and the selection of primaries and satellites in § El In §m 
we derive the interloper-corrected surface density pro¬ 
file of satellites in volume-limited and flux-limited SDSS 
samples, and their subsamples, and compare results to 
the ACDM cosmological simulations. We also discuss 
comparisons to simulations results. Our main results 
and conclusions are summarized in section Through¬ 
out this paper, we assume flat ACDM cosmology with 
fl m = 0.3 and h = 0.7. 

2. INTERLOPER SUBTRACTION 

We use cosmological A'-body simulations to construct 
mock samples of host halos and possible satellites, with 
which we can examine the effects of interloper contam¬ 
ination on the satellite distributions and test different 
methods of interloper subtraction. We start with two 
80ft -1 Mpc high-resolution dark matter simulations of 
the concordance ACDM cosmology: H m =0.3, h= 0.7, and 
(78=0-9- The two simulations differ in the random seed 
of their initial conditions. We use outputs of the simu¬ 
lations at a redshift of 0.1, which is similar to the aver¬ 
age redshift of objects in the SDSS spectroscopic survey. 
The simulations were performed w ith the A da ptive Re¬ 
fin ement Tree (ART) IV-body code (iKravtsov et al.ll997i: 
|Kr avtsojll999|). Details o f the simulations can be found 
in ITasitsiomi et ail (12004 . Halo identification was per¬ 
formed using a variant of the Bound Dens ity maxima 
halo finding algorithm llKlvnin et alJll999all . Details of 
the algorith m and parameters used in the halo finder can 
be found in IKravtsov et all ( 2004al ). 

In addition, we use a simple, observationally motivated 
scheme to assign luminosi ties to the halos. Details of this 
method can be found in ITasitsiomi et a.U 1 2004 1. Ha¬ 
los are assigned luminosities by matching the cumulative 
velocity function, n(> Umax), where U max is the maxi¬ 
mum circular velocity, to the SDSS luminosity function, 
n{< M r ) at z = 0.1 ([Blanton et a,l J l2003rf) . M r is the 
SDSS r-band absolute magnitude defined as M — 51og/i. 
The magnitudes have been K-correcte d to z=0.0, using 
kcorrect v3.2 ijRlanton et a,l .ll2003rfl . Scatter is intro¬ 
duced in the relation between Unax and M r , assuming 
a standard deviation of 1.5 magnitudes for the M r dis¬ 
tribution at fixed U max - All galaxies down to M r = -18 
(corresponding to a mean minimum Unax = 100 km s -1 ) 
are included. 

In the remainder of the paper, our terminology is to re- 
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TABLE 1 

Selection & Isolation Criteria for Test Samples 


Parameters 

Test Sample 1 (TS1) 

Test Sample 2 (TS2) Test Sample 3 (TS3) 

Constraints on primaries 

Vmax = 100-150, ..., 300-350 km s” 1 

Mr < -20 

Mr < -20 

Satellite objects 

DM particles 

DM particles 

satellite galaxies 

Isolation criteria: 

Size difference 

Wax > 0.5Vmax 

Wax > 0.5Vmax 

A Mr < 2 

Minimum projected distance, AR(h ~ 1 Mpc) 

0.5 

0.5 

0.5 

Minimum velocity separation, AV (km s —1 ) 

1000 

1000 

1000 

Satellite sample criteria: 

Magnitude difference from the primary 

— 

— 

A Mr > 2 

Maximum projected distance, Sr(h ~ 1 Mpc) 

0.6 

0.6 

0.6 

Maximum velocity separation, Sv (km s —1 ) 

500 

500 

500, 1000 

Number of isolated primaries 

380, 289, 236, 143, 89 

728 

728 

Number in satellite sample 

7475, 11165, 16614,15354, 14899 

50608 

343, 401 

Limiting magnitude M r 

— 

— 

-18 


fer to the sample of hosts or primaries constructed using 
isolated DM halos (i.e., objects which do not lie within 
a virial radius of a larger object) selected by their max¬ 
imum circular velocities as a sample of “halos” and the 
sample of isolated DM halos assigned r-band luminosi¬ 
ties as a sample of “galaxies,” respectively. We use V max 
to quantify the size of halos because it is measured more 
robustly and not subject to the same ambiguity as mass 
definitions. For “galaxies”, we use magnitudes, M r , to 
quantify the size of galaxies, as an alternative to V max to 
account for possible effects of scatter between V max and 
M r . 

For possible satellites, subhalos assigned r-band lu¬ 
minosities are referred to as ’satellite galaxies’ and are 
selected by their M r magnitude difference from their 
parent ’galaxies.’ For isolated halos, satellites are se¬ 
lected from an arbitrary fraction of dark matter particles 
around each host halo. These two choices of test satel¬ 
lite samples are expected to bracke t the possib l e range of 
radial profiles of the rea l satellites llDiemand et alJl2004 
iNagai fc Kravtso\l20f) 5A Below we detail the definition 
of samples in our analysis. 

2.1. Primary and Satellite Samples 

Previous studies of the projected distribution of galac¬ 
tic satellites used observationally motivated selection cri¬ 
teria to construct a primary sample of isolated host 
galaxies and a sample of potential satellites that are 
projected close to primaries. To test several alternative 
methods for interloper subtraction, we construct a set 
of primary and satellite samples derived from simula¬ 
tions using different isolation criteria for the primaries 
and different selection of the satellite samples. 

For th e primaries we use isolation criteria similar to 
those of iPrada et al . (2 1)l)3i . We start by creating pri¬ 
mary samples of halos and galaxies. An isolated primary 
halo, with circular velocity yP” x , must have no other ha¬ 
los with a maximum circular velocity Vmax > 0.5 V^ax 
(which corresponds to an absolute magnitude difference 
of ss 1.5 — 2 magnitudes) within a projected separa¬ 


tion A R = 0.5 h~ x Mpc and velocity separation, AV = 
1000 kms -1 . An isolated galaxy primary with absolute 
magnitude Mj? rl must have no other galaxies brighter 
than MP rl + 2, within the same projected separation and 
velocity separation. 

For each isolated primary, we construct two satellite 
samples. In the first, we use a random subset of DM 
particles found within a projected distance, Sr = 0.6 
ft , -1 Mpc, and velocity difference, Av = 500 kms -1 , 
from each primary. These are our fiducial choices of Sr 
and Sv. We also test a non-fiducial velocity difference, 
Sv = 1000 kms -1 . In the second satellite sample, we 
select all galaxies fainter than the primary by more than 
two magnitudes, within the same projected distances and 
velocity difference as above. 

For our primary halos, we create several DM particle 
satellite samples for different ranges of VjP” x (100-150, 
150-200, 200-250, 250-300, 300-350 km s -1 ). This set of 
DM halo primaries and DM particle satellites is referred 
to hereafter as Test Sample 1. For primary galaxies, we 
use a single sample of 728 galaxies with M r < — 20 and 
build three different satellite samples using DM particles 
with the Sv = 500 km s -1 criterion, satellite galaxies 
with Sv = 500 km s -1 , and satellite galaxies with Sv = 
1000 km s -1 . The set of primary galaxies and DM parti¬ 
cle satellites is referred to as Test Sample 2, while the set 
using satellite galaxies - subhalos with assigned luminosi¬ 
ties - and both velocity criteria is labeled Test Sample 3. 
These samples are summarized in Tabled For samples 
with DM satellite particles, we bin the objects in radial 
bins of 35/i -1 kpc, starting at a minimum separation of 
25.6/i -1 kpc. For samples with satellite galaxies, we bin 
objects in bins of 70/i -1 kpc, starting at the same mini¬ 
mum separation. The larger bin in the latter case is due 
to the smaller statistics of the subhalo satellite sample. 

2.2. True Satellites vs. Interlopers 

There is a fraction of objects in our satellite samples 
that are not gravitationally bound to the primaries but 
are included in the sample because of projection effects. 
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Fig. 1.— Testing interloper subtraction with mock catalogs: 
the derived average projected distribution of satellites compared to 
true satellites for the mock samples with five V^ax ranges: 100-150, 
150-200, 200-250, 250-300, and 300-350 km s -1 (Test Sample 1). 
Top: The projected radial number density profile for true satellites 
(solid) and the satellite sample (dotted). The amplitude in number 
density increases with V5Jax- Center: The fraction of interlopers 
in the satellite sample as a function of projected radius. The bins 
are distinguished by line type: 100-150 (solid), 150-200 (dotted), 
200-250 (short dashes), 250-300 (long dashes), and 300-350 (dot- 
dashed) km s -1 . Bottom: The projected interloper number density 
profile using the same line types as the center panel. 


Throughout this paper, we call such objects interlopers. 
We define true satellites as objects which satisfy the fol¬ 
lowing negative binding energy criterion: 

E < ^|F| 2 + <K r )> where (j){r) = — Vesc ^ . (1) 

Here V = v sat — v pr i is the 3D satellite velocity with 
respect to the primary, and, for a iNavarro et al . con 
hereafter NFW) profile, 


v 2 esc (r) = 9.26 V; 


2 ln(l + r/r s ) 


r/r s 


( 2 ) 


where r s has been estimated as r a = r max /2.15 and r max 
is the radius of the maximum circular velocity of the host 
halo, Vmax, as expected for a NFW density profile. 

Using our test samples, we look for biases in the true 
satellite surface density profile and the projected distri¬ 
bution of the satellite sample. We calculate the surface 
density of satellites per primary by 


(E(i?)} = 


-Win 

nN pii (Rl-RlY 


(3) 


where Win is the number of objects in the satellite sam¬ 
ple that are found between the inner, i?i and outer, R?, 
radii of the annulus and 7V P ri is the number of objects in 
the primary sample. R is the midpoint of the bin. 

The projected number density profiles of all satellites 
and only true satellites as a function of distance to the 
primary for Test Sample 1 are plotted in the top panel of 
Figure n Here we can see that interlopers significantly 
flatten the projected radial distribution. The flattening 
is the strongest for the surface density profiles of satellite 
samples of the smallest primary halos. For the different 
UP” primary samples, the DM particle satellite samples 
show that the fraction of interlopers in the satellite sam¬ 
ple as a function of projected radius is similar for all mass 
ranges, but decreases with increasing primary mass. In 
addition, at all masses, the projected number density of 
interlopers (bottom panel of Fig. |T|| is relatively flat, but 
rises at radii greater than the isolation criterion. The re¬ 
sults for Test Samples 2 and 3 are similar to those shown 
in Figure Q 


2.3. Interloper Subtraction Methods 
2.3.1. Random Points Method 

The most straightforward procedure for interloper re¬ 
moval, which we call the random points method, assumes 
that the surface density of interlopers is constant and 
can be estimated b y surveying the area arou nd randomly 
placed points Ce.g.. lLake fc TremainellTflSOl) . This is one 
of the simplest methods to apply, since it can be imple¬ 
mented even without redshift information. 

We apply this method to our simulation-derived sam¬ 
ples, by choosing random points within the simulation 
box which satisfy the same isolation criteria as the corre¬ 
sponding primary sample. Each of these random points, 
which make up a mock primary sample, has the charac¬ 
teristics - mass, magnitude, and velocity - of an object 
in the real sample of primaries. This is accomplished by 
sampling the primary sample with replacement - i.e., as¬ 
signing each random point with the characteristics of an 
object in the primary sample, choosing that object ran¬ 
domly from the entire primary sample. To estimate the 
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Fig. 2.— The fraction of primaries with a number of objects 
with Vjnax > 0.5Vmax hi the annulus of 0.5 < A R <1.0/i _1 Mpc 
and Ay < 500 kms -1 . The primary sample of y m ax = 200-250 
km s _1 (Test Sample 1) is plotted here (thin, solid line) and the 
corresponding random points mock primary sample is shown in the 
thick, dashed line. 



AV (km s -1 ) 


Fig. 3.— Testing interloper substructure with mock catalogs: the 
relative velocity distribution of the interloper sample in the mock 
catalog for primaries with V m ax = 200 - 250 km s _ 1 (Test Sample 
1) for all DM particles within 25.6 < A R < 500/i —1 kpc (solid line, 
Poisson errors). The estimated interloper sample for the clustered 
random points method is shown as a dotted line and for the random 
points method as a dashed line. The nearby points methods are 
shown in dot-dashed lines, with square points for Ait corr = 1 - 
2h ~ 1 Mpc and crosses for Ai? corr = 1 - 5/v 1 Mpc. The thin 
vertical line shows the velocity difference criterion of the satellite 
sample of 500 km s -1 . 


interloper contribution, we construct a sample of random 


points that has 20 times the number of objects in the 
sample of primaries. We can expect, however, that uni¬ 
formly distributed random points will preferentially sam¬ 
ple voids rather than the environments similar to those of 
the isolated primaries. This is because 1) voids occupy 
most of the volume and 2) our isolation criteria have 
forced the random points to preferentially lie in voids, 
because the fraction of space that satisfies our isolation 
criteria is larger in voids than in dense environments. 

This is illustrated in Figure [21 which shows the distri¬ 
bution of massive (V max > 0.5V^x) neighbors within 
the annulus described by the projected separation of 
0.5 < A R < 1.0 h~ l Mpc (i.e., just outside the radius 
used in the isolation criteria) and velocity difference of 
AV = 500 kms -1 for our sample of primaries (Test 
Sample 1) and the sample of random points satisfying 
the same isolation criteria. The Figure [21 shows that 
primaries have systematically larger number of massive 
neighbors than the random points, and, therefore, that 
primaries occupy clustered regions of space. The random 
points method is thus expected to consistently underes¬ 
timate the interloper contamination. This result may 
also be surmised from the bottom panel of Fig. Q which 
shows that the number density of interlopers is not uni¬ 
form and is correlated with primary halo mass. More 
massive isolated halos have more interlopers, which is 
consistent with the concept that, in general, more mas¬ 
sive objects are more clustered. 

2 . 3 . 2 . Clustered Random Points Method 

It is clear that we can improve upon the random points 
method, if we select not random points but points that 
are in environments similar to those of primary galax¬ 
ies. We attempt to achieve this in the clustered random 
points method. Isolated points are chosen at random and 
the number of objects in the annulus described above is 
calculated. Points that allow the fractional distribution 
of annular objects from the clustered random points to 
exceed that of the primary points are rejected. In other 
words, the method insures that the mock sample con¬ 
tains the same distribution of massive neighbors around 
the primaries, as shown in Figure [21 As we discuss be¬ 
low, the clustered random points method performs con¬ 
sistently better than the random points method in esti¬ 
mating the interloper contamination. 

In Figure [21 we show the relative velocity distribu¬ 
tion of interlopers and estimated interlopers for primaries 
with Vjnax = 200 — 250 kms -1 in Test Sample 1. The 
interloper distribution is not constant as a function of 
velocity, showing a peak at AV = 0 and a tail to AV 
= 1000 km s -1 . The interloper population includes not 
just objects uniformly distributed in the velocity space, 
but also objects with velocities correlated with the veloc¬ 
ity of the primaries. The estimated interloper population 
using the random points is uniformly distributed in veloc¬ 
ity space and significantly underestimates the interloper 
contribution (dashed line). The clustered random points 
method fares considerably better, although the number 
of interlopers is still somewhat underestimated (dotted 
line). 

2 . 3 . 3 . Nearby Points Method 

An alternative way of ensuring that the random points 
sample environments of the primaries correctly is to pick 
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isolated points that are within a pro jected correlation 
length of real primaries. (Smith et al . I2004t l use such a 
method, estimating the background using points at pro¬ 
jected distances > 350/i -1 kpc. In their isolation crite¬ 
ria, they require that the magnitude difference between a 
neighbor and the primary must be greater than 0.7 mag¬ 
nitudes for galaxies within a projected distance of 700 
h~ x kpc. The area outside of a 350/i -1 kpc projected ra¬ 
dius, however, may not have the same isolation criteria 
as the primary galaxy sample. 

Observational measurements of the two-point correla¬ 
tion function of bright galaxies find a co rrelation length 
of 5 h^ 1 Mpc fe.g.. iZehavi et all 120041. In Figure 0 
we show the estimated fraction of interlopers using a 
method in which test points are selected from the an¬ 
nuli of A R COII = 1 — 2h~ 1 Mpc and AR corr = 1 — 5 h _1 
Mpc around the primaries. We choose the inner radius 
of 1 h~ x Mpc to avoid sampling real satellites. All test 
points satisfy the same isolation criteria as our sample 
of primaries. The Figure 0 shows that the velocity dis¬ 
tribution of interlopers for AR corr = 1 — 5 h _1 Mpc is 
somewhat similar to that of the clustered random points 
method, while AI? corr = 1 — 2 h~ x Mpc choice recov¬ 
ers the true interloper velocity distribution much better. 
We will, therefore, use this latter radial annulus as our 
fiducial choice. As we show below, the nearby points 
method with the fiducial A R co „ is the best among the 
other methods we tested here in recovering the interloper 
contamination in our mock samples. 

2.3.4. ’’Gaussian + Constant” Method 



AV (km s -1 ) 

Fig. 4. — The Gaussian plus constant method (thick, solid line) is 
compared to the velocity distribution of the satellite sample (square 
points, dotted line) and the interloper sample (dashed line) for 
primaries with 14iax = 200-250 km s —1 and satellites with 25.6 < 
A R < 500h —1 kpc (Test Sample 1). The thin horizontal line shows 
the estimated interloper fraction from the Gaussian plus constant 
method. 

For completeness, we also test th e interloper sub- 
traction method used by IMcKav et all {2002) and 
iPrada et all (I2003fl . The method does not use random 
points, but assumes instead that the velocity distribu¬ 
tion of satellites can be described by a Gaussian, while 


the distribution of interlopers is uniform. The velocity 
distribution of interlopers shown in Figure 0 is inconsis¬ 
tent with this assu mption; a resu l t that i s in a greement 
with conclusions bv lvan den Bosch et ~an i(200l see their 
Fig. 1). 

Nevertheless, the tests performed on our mock samples 
show that the method does estimate the velocity disper¬ 
sion of the true satellites - t he purpose for which the 
method was origin ally used hv lMcKav et all (j20D2ll and 
IPrada et all 1200311 - quite accurately. This is because 
the velocity dispersion of the interlopers, correlated with 
primaries in velocity space, is similar to that of the satel¬ 
lites. Their inclusion into the satellite samples thus does 
not bias the velocity dispersion appreciably. The veloc¬ 
ity dispersion of interlopers is likely due to the infall of 
objects along filaments. 

Nevertheless, to correctly account for the interloper 
contribution to the radial surface density profile this is 
not sufficient. The assumption that all of the objects cor¬ 
related with the primary in velocity space are true satel¬ 
lites will lead to an underestimate of interloper contribu¬ 
tion. This is demonstrated in Figure Q] Given that the 
interloper fraction increases with increasing projected ra¬ 
dius, the ’’Gaussian plus constant” method leads to sur¬ 
face density profiles flatter than the true distribution. 
Another problem occurs when the number of satellites is 
small and their velocity distribution is not well sampled. 
We find that in such cases the Gaussian+constant fits 
can be unstable. 

In our tests, we use objects within Sv = 1000 kms -1 
to fit a gaussian+ constant. This is larger than Sv = 
500 kms -1 used bv lPrada et al . ( |2003ll and in this paper 
fiducially, because our sample includes massive objects 
with velocity dispersion close to the Sv, in which case 
the interloper constant is poorly constrained. 

2.4. Testing the Interloper Subtraction Methods 

We compare how well the four different interloper sub¬ 
traction methods described above recover the true pro¬ 
jected density profile of satellites (E(i?)) true , where 

<E(i2)> int sub = (£(£)> sat - <£(-R))est int, (4) 

where our interloper subtracted density profile, 
(£(!?)}i n t sub, is the profile of the satellite sample, 
(£(!?)) sat, minus the interloper profile, (£(i?.)) es t int, 
estimated using different methods described above. 

Figure 0 shows the results for Test Sample 1, for dif¬ 
ferent ranges of our primaries. The figure shows 
that the random points method performs the worst. The 
’’Gaussian plus constant” and clustered random points 
methods perform better, although the latter does a bet¬ 
ter job overall in estimating the interloper contamina¬ 
tion. The nearby points performs best. For Test Sample 
2 and 3, the results are similar, except for the ’’Gaus¬ 
sian plus constant” method, which runs into difficulties 
with fitting correct interloper level for Test Sample 3, in 
which the number of satellites is small and the velocity 
distribution is under-sampled. This indicates that the 
method should be applied only to the samples in which 
the velocity distribution of satellites and interlopers is 
sampled well. 

2.5. Fits to the radial profiles 
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Fig. 5.— Testing interloper subtraction with mock catalogs: the 
fraction of interlopers for different mock satellite samples estimated 
by different methods. Clockwise from the top, left panel, the Vlnax 
ranges are 150-200, 200-250, 250-300, and 300-350 km s —1 , using 
DM particle satellites (Test Sample 1). The corresponding number 
density profiles are shown in Fig. Q(top panel). In all plots, the 
thick lines represent the true interloper fraction (also shown in Fig. 
HI center panel), the dotted lines use the Gaussian plus constant 
method, the dashed lines use the random points method, the solid 
lines use the clustered random points method, and the dot-dashed 
lines use the nearby points method, with Aitcorr = 1 — 2h~ 1 Mpc. 



We now test how well the fits to the interloper sub¬ 
tracted radial surface density profiles recover the profile 
of the true satellites. The data is binned with bins of 
35 h~ l kpc starting at a projected radius of 25.6 h -1 kpc 
for the DM particles as satellites, and with bins of 70ft -1 
kpc for subhalos with assigned luminosities starting at 
a projected radius of 25.6h -1 kpc. We fit the number 
density profile with a power law, 

T,{R)=AR a , (5) 

for bins with mean projected radius smaller than 0.5/i _1 
Mpc. We find the best fit values of the slope, a, and its 
errors marginalized over the normalization of the power 
law, for the true satellite sample, the satellite sample 
(true satellites + interlopers), and for the interloper sub¬ 
tracted profile obtained with each of the methods dis¬ 
cussed above. We measure the ler errors in slope, which 
reflects the statistical Poisson errors of the density pro¬ 
files. 

In general, the true DM particle satellite distribution 
should be well-approximated by a projected NFW pro¬ 
file. However, we do not attempt to fit an NFW profile 
because, as we can see below, the current observational 
samples cannot discriminate between the NFW and a 
simpler power-law profile. We choose our sample of pri¬ 
maries (Test Samples 2 & 3) using luminosities instead 
of Vjnax and stack the satellites of many primaries, pos¬ 
sibly mixing NFW profiles with different concentrations 
and virial radii together. In addition, due to the limited 
satellite statistics and the fiber collisions in the SDSS we 
use a minimum projected radius, similar to that used in 
the SDSS spectroscopic sample (see below), which is of 


order of 10% of the virial radius of the primary galaxies. 
The minimum radius, then, may be larger than the scale 
radius of the NFW profile. 

As can be seen from Tables [S] and |3j the DM particle 
satellite sample (Test Sample 2) and the galaxy satellite 
sample (Test Sample 3) have different radial distribu¬ 
tions. Most of the DM particles in the satellite sample 
are from the smooth distribution of the parent halo, not 
subhalos, which explains the steep, a = —1.81 ± 0.01, 
slope found for the DM particle satellite population. This 
is consistent with the steep slope of ~ —3 (or ~ —2 in 
projection) predicted for the 3-D density profile of the 
CDM halos at large radii. The radial distribution of the 
satellite sample which uses subhalos with assigned lumi¬ 
nosities is considerably flatter a = —1.34±0.12. At small 
projected radii, the profile is flattened by tidal disruption 
of subhalos. 


TABLE 2 

Estimated Power-Law Slope for Sample with DM Particles 
(Test Sample 2) 


Input Data 


25.6 < R < 500 h- 1 kpc 


true satellites 
satellite sample 
random points 
clustered random points 
nearby points (ATJcorr = 
gaussian + constant 


-1.815 ±0.007 
-1.352 ±0.007 
-1.481 ± 0.008 
-1.570 ±0.008 
1 - 2) -1.686 ±0.009 

-1.533 ±0.007 


TABLE 3 

Estimated Power-Law Slope for Sample with Satellite 
Galaxies (Test Sample 3) 


Input Data 


25.6 < R < 500 h 1 kpc 


8v = 500 km s 1 


true satellites 


-1.34 

± 

0. 

12 

satellite sample 


-0.84 

± 

0. 

11 

random points 


-0.95 

± 

0. 

12 

clustered random points 


-1.08 

± 

0. 

12 

nearby points (AR corr = 1 - 

2) 

-1.20 

± 

0. 

14 

gaussian + constant 


-0.95 

± 

0. 

11 

' = 1000 km s -1 






true satellites 


-1.34 

± 

0. 

12 

satellite sample 


-0.75 

± 

0. 

11 

random points 


-0.95 

± 

0. 

12 

clustered random points 


-1.23 

± 

0. 

14 

nearby points (AR CO tt = 1 - 

2) 

-1.24 

± 

0. 

15 

gaussian + constant 


-0.94 

± 

0. 
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In Tables 0 and 0 we show that the best-fit slopes 
for the satellite radial profile are significantly flattened 
by interlopers. Without any interloper subtraction, the 
estimated surface density profile will be shallower than 
the profile of the true satellite population, in both the 

















DM particle sample and the galaxy sample by ss 0.5 in 
the power-law slope. 

The random points method is inadequate for recov¬ 
ering the correct slope: the recovered slope slopes are 
shallower than those of the true distribution by « 0.4. 
The ’’Gaussian plus constant” method also underesti¬ 
mates the slope significantly. Interloper subtraction by 
the clustered random points and nearby points perform 
reasonably well in all test samples. In general, the 
nearby points method performs the best giving on aver¬ 
age a steeper slope (by ~ 0.1) than the clustered random 
points method. Overall, the interloper bias still persists 
as the slope estimated with the nearby points method 
systematically underestimates the true slope by ~ 0.1, 
even though the slope values estimated with this method 
are within one standard deviation of the true slope for 
Test Sample 3. We will therefore use the nearby points 
method as our interloper correction method of choice, 
keeping in mind that the best fit slope should be cor¬ 
rected by Aabias ~ 0.1. 

3. THE SDSS SPECTROSCOPIC SURVEY 

The Sloan Digital Sky Survey (SDSS) ( [York et al . 
Hoi will image up to 10 4 deg 2 of the northern 
Gala ctic cap in five bands, u,q,r,i.z, down to r ~ 

22.5 llFukugita et al .11 1 9961iHogg et aLll2001l: iSmith et al.l 

2002) using a dedicated 2.5 m telescope at Apache P oint 
Observatory in New Mexico (IGnnn et al .111 Eio3l . In 
addition to the imaging survey, the SDSS main galaxy 
sample is a subsample of objects from the imaging cata¬ 
log w hich have been targeted for spectroscopic observa¬ 
tions llStrauss et al.ll2f)02D . The spectroscopic targets are 
selected with r-band Petrosian magnitudes r < 17.77 and 
r-band Petrosian half-light surface brightnesses < 

24.5 mag arcsec -2 . The median redshift of the SDSS 
main galaxy sample is 0.104. 

The SDSS spectroscopy is carried out using 640 optical 
fibers positioned in pre-drilled holes on a circular plate in 
the focal plane of 3 degrees diameter, with minimum sep¬ 
aration between fibers of 55". Targeted imaging regions 
are assig ned spectroscop ic plates by an adaptive tiling al¬ 
gorithm llBlanton et alJl2003b ). which also assigns each 
object a fiber. An automated pipeline meas u res t he red- 
shifts and classifies the reduced snectra llStoughton et al.l 

f200a; pier et al MM llvezic et al.12004 D. J. Schlegel et 

al. 2005, in preparation). 

For this catalog we use the reductions of the SDSS 
spectroscopic data performed by D. J. Schlegel et al. 
(2005, in preparation) using their reduction code, which 
extracts the spectra and finds the redshifts. The red- 
shifts found are, over 99% of the time for Main galaxy 
sample targets, identical to the redshifts found by an al¬ 
ternative pipeline used for the SDSS Archive Servers (M. 
SubbaRao et al. 2005, in preparation). 

For this analysis, we use a subset of the available spec¬ 
troscopic main galaxy sample released as of Data Re¬ 
lease Four ( Adelman-McCarthv et akfeOO B!. but includ¬ 
ing all of the galaxy sample r eleased as of Data Release 
Three (lAha.zaiia.n et a,l.|l2005 1. This catalog, known as 
LSS SAMPLE14, is built from th e New York Univer- 
sity Value-Added Galaxy Catalog l|Blanton et al .1 [20051) 
and contains 312,777 galaxies. Because the SDSS spec¬ 
troscopy is taken through circular plates with a finite 
number of fibers of finite angular size, the spectroscopic 


completeness varies across the survey area. The result¬ 
ing spectroscopic mask is represe nted by a combination 
of disks and spherical polygons llTegmark et al.l [20041. 
Each polygon also contains the completeness, a number 
between 0 and 1 based on the fraction of targeted galax¬ 
ies in that region which were observed. We apply this 
mask to the spectroscopy and include only galaxies from 
regions where the completeness is greater than 90%, an 
area of 3448 square degrees. The same criterion is ap¬ 
plied for catalogs of clustered random points and nearby 
points used for the interloper subtraction. 

4. RADIAL DISTRIBUTION OF SATELLITES AROUND 
SDSS PRIMARIES 

Unlike the numerical simulations where we have good 
resolution and 100% completeness, when using spectro¬ 
scopic data to stack objects and estimate a surface den¬ 
sity, we suffer from two major problems: fiber collisions 
and incompleteness. 

As described previously, the minimum separation be¬ 
tween fibers, the fiber collision separation, is 55". At 
a redshift of z = 0.035, the approximate median of our 
SDSS sample of primary galaxies, the fiber collision sep¬ 
aration is 26.8 ft -1 kpc. Fiber collisions could bias the 
small projected radii end of the radial distribution, re¬ 
moving objects that should be counted and tilting it shal¬ 
lower. However, some of the area of the survey has been 
re-observed and the overlap region could have objects 
observed with separations as small as the fiber diameter, 
3". In our samples, we do not use objects at projected 
radii smaller than the fiber collision separation. 

We use the r-band magnitudes in the LSS SAMPLE14 
subsample of the SDSS spectroscopic survey, normalized 
to h= 1, such that M r = M 0 .i r — 51og 10 /i, where M 0 .i r 
is the abs olute magnitude K-co rrected to z—0.1 as de¬ 
scribed in lBlanton et al.1 ( 2003a !. The LSS SAMPLE14 
also provides measures of the fraction of objects with 
spectra in the area of an object. While fiber collisions 
remove objects from our survey, incompleteness removes 
area from our survey, which would likely tilt the number 
density steeper, since there is more area further away 
from a primary galaxy. Constraining our sample to in¬ 
clude only galaxies with a minimum completeness frac¬ 
tion of 90% should lead us to be complete to that level. 
However, since we need to search the area around each 
primary galaxy, it is possible that the search area will not 
be contained on an area of one completeness level but can 
overlap with an area of a lower completeness percentage. 
A simple check of this is to calculate an analogue to the 
projected cross-correlation function, w(R ), which is un¬ 
biased by incompleteness and can be compared to the 
projected number density estimated using the random 
points method. We apply such test and show that in¬ 
completeness does not bias our estimates of the radial 
profiles (see §E3- 

We create both a vol ume- limited sample (§ ITU and 
a flux-limited sample (§ Oil . The flux-limited sample 
offers better statistics, but it is also biased to brighter 
satellite objects and the radial profile of satellite galax¬ 
ies has not been established to be independent of the 
magnitude of satellites. For this reason, we also test a 
volume-limited sample with poorer statistics. 

4.1. Volume-limited, samples 
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TABLE 4 

Selection & Isolation Criteria for SDSS Samples 


Parameters 

Volume-Limited 

Flux-Limited 

Maximum depth of sample (km s —1 ). 

13,500 

13,500 

Constraints on primaries. 

-23 < M r < -20 

-23 < M r < -20 

Constraints on bright neighbors: 

Maximum magnitude difference. 

2 

2 

Minimum projected distance, AR(h~ 1 Mpc).... 

0.5 

0.5 

Minimum velocity separation, AV (km s —1 ). 

1000 

1000 

Constraints on satellites: 

Minimum magnitude difference. 

2 

2 

Maximum projected distance, Sr(h —1 Mpc). 

0.5 

0.5 

Maximum velocity separation, Sv (km s — - 1 ). 

500, 1000 

500, 1000 

Minimum projected distance, Sr{h~ 1 Mpc). 

0.329, 0.341 

0.329, 0.341 

Number of isolated primaries. 

871 

871 

Number in satellite sample. 

336, 357 

678, 786 

Limiting magnitude M r . 

-17.77, -17.85 

— 


300 

| 200 
1 

Z 100 

0 




M 


Satellites per primary 
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Fig. 6. — Statistics of primaries and satellites for a volume- 
limited sample. Top-left: The r-band magnitude histogram for 
primaries. Bottom-left: The r-band magnitude histogram for satel¬ 
lites. The solid line shows the results for the Sv = 500 km s —1 
criterion, while the dotted line shows the Sv = 1000 km s —1 crite¬ 
rion. Bottom-right: The magnitude differences between satellites 
and primaries, with line styles as in the bottom-left panel. Top- 
right: The number of satellites per primary for primaries with at 
least one satellite, with line styles as in the bottom-left panel. 


We first create a volume-limited galaxy sample from 
the LSS SAMPLE14 with a depth of 13,500 km s -1 , cor¬ 
responding to the limiting redshift of z = 0.045. This 
limit is chosen as a trade-off between the volume of the 
sample and the absolute magnitude limit for our satel¬ 
lites, which would need to be decreased to brighter mag¬ 
nitudes for more distant primaries. The trade-off is also 
with the minimum separation at which fiber-collisions 
become important, which increases with distance. To in¬ 
clude more distant primaries we would have to sacrifice 


the ability to probe density distributions at small sepa¬ 
rations. 

In total we have 21,851 galaxies. Since the isola¬ 
tion criterion requires that we search for objects that 
are within AV = 1000 km s -1 , we can only search 
for primaries within the subset of velocities 1000 to 
12,500 km s -1 . For the satellite catalog to be volume- 
limited, this requires a maximum absolute magnitude of 
A4r,iim — 51og/i = 17.77 — DM — Kq.i in the r-band, where 
the 17.77 is the flux limit in this band, DM is the dis¬ 
tance modulus, and Ao.i is the K-correction at 2=0.1. 
We use the A'-correction at 2 = 0.1 in order to avoid un¬ 
derestimating the limiting absolute magnitude. As in the 
simulations, we test both the Sv = 500 km s -1 satellite 
criterion and a larger Sv = 1000 km s -1 . For those limits, 
the minimum separation between fibers is 32.9 h~ l kpc 
and 34.1 h _1 kpc, respectively and the limiting absolute 
magnitudes are -17.77 and -17.85. The satellites are thus 
limited to the brightest satellite galaxies, ~ 0.1A*. We 
choose galaxies that are in areas that are at least 90% 
complete and set the size of the mock primary sample to 
be 20 times the number of primary galaxies. 


TABLE 5 

Estimated Power-Law Slopes for the Volume-Limited 
Sample 


Input Data 

Sv = 500 km s —1 





satellite sample 


-1.21 

± 

0.09 

clustered random points 


-1.46 

± 

0.11 

nearby points (AR CO rr = 1 - 

2) 

-1.58 

± 

0.11 

Sv = 1000 km s —1 





satellite sample 


-1.18 

± 

0.09 

clustered random points 


-1.55 

± 

0.11 

nearby points (AR CO rr = 1 - 

2) 

-1.65 

± 

0.12 


The statistics of primaries and satellites for both sam¬ 
ples are shown in Figure|SJ The number of possible satel¬ 
lites found in the volume-limited samples is small. For 
the range —23 < M r < —20, there are 871 primary galax- 
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Fig. 7.— The satellite sample compared to interloper subtracted 
samples for the volume-limited sample using 8v = 500 km s —1 . 
Top: The best-fit power-law for the satellite sample (solid) and 
the interloper subtracted profile for the AR CO rr = 1-2 h~ 1 Mpc 
nearby points (dotted). The data is shown with error bars: satellite 
sample (square points, thick Poisson error bars) and nearby points 
(open circles, thin error bars). Bottom: The 68%, 90%, and 99% 
confidence intervals for the satellite sample (solid lines) and the 
AR COII = 1-2 h- 1 Mpc nearby points sample (dotted lines). The 
contours are labeled by the unmarginalized best-fit slope of each 
profile. 


ies and 336 objects in the satellite sample with projected 
radii greater than the minimum separation and less than 
0.5 ft -1 Mpc. For the Sv = 1000 km s -1 sample, which 
uses a slightly larger minimum separation and limiting 
magnitude, there are 357 galaxies. The volume-limited 
samples are summarized in Tabled 
The satellite sample and the nearby points interloper 
subtracted results are shown in Figure 0 for the veloc¬ 
ity criterion of Sv = 500 km s -1 , in bins of 70 h^ 1 kpc, 
starting from the minimum separation of 32.9 h^ 1 kpc. 
The results are similar to those in the simulations; the 
satellite sample distribution is shallower than the inter¬ 


loper subtracted samples. In addition, the nearby points 
method distribution is steeper than that of the clustered 
random points method. 

We fit the radial profile with a power law, with the 
results for the slope marginalized over the amplitude of 
the power law shown in Table 0 The slope of the best- 
fit power law for the satellite sample is —1.21 ± 0.09 in 
the sample with Sv = 500 km s^ 1 , with very similar 
results for the Sv = 1000 km s -1 sample. The clustered 
random points method finds a slope of —1.46 ±0.11 and 
the nearby points methods finds a slope of —1.58 ± 0.11 
in the Sv = 500 km s _1 sample. The Sv = 1000 km s _1 
sample shows steeper fits, with slopes of —1.55 ± 0.11 
and 1.65 ± 0.12, respectively. The systematic differences 
between the best-fit slopes of the clustered random points 
and nearby points is ~ 0.1, consistent with the result 
found in the simulations. In the simulations, the bias 
found in the interloper subtraction methods was ~ 0.1 
for the nearby points method and ~ 0.2 for the clustered 
random points method, which would imply that the slope 
of the true satellite distribution is a = a e st — Actbias ~ 
-1.7. 

The marginalized errors suggest that slopes of the 
satellite samples and the interloper subtracted samples 
are significantly different. We illustrate this point in Fig. 
0 where we plot all the points at projected radii smaller 
than 0.5 h^ 1 Mpc for the Sv = 500 km s , satellite sam¬ 
ple and the results of the nearby points method. Here the 
best-fit power-law is plotted in the top panel, where the 
best-fit slope of the satellite sample is a = -1.22 and the 
best-fit slope of the interloper subtracted sample is a =- 
1.58. The bottom panel shows the confidence regions for 
the two fits, where the slopes of the two distributions do 
not overlap within the 99% confidence intervals. 

4.2. Completeness test 



Fig. 8.— The ratio of the cross-correlation function, w(R), to the 
projected, interloper subtracted surface density of true satellites 
estimated using the random points method for the volume-limited 
sample with 5v = 500 km s —1 . 

Although we have selected our sample carefully to 
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avoid spurious results, we still need to assess to what 
extent our results might be biased by sample incom¬ 
pleteness, caused by, for example, holes in the survey 
or edge effects. We can test this by calculating for our 
primary-)- satellite sample a statistic analogous to the 
projected cross-correlation function: 


w(R) = 


(N sa , t (R)) 

(N int (R)) 


( 6 ) 


where {N sat (R)) = (IVtruesat(-R)) + (Mnt {R)) and 
{Ni n t(R)) is the estimated number of interlopers at sep¬ 
aration R in a sample, and (N truesa t(R)) is the corre¬ 
sponding estimated average projected number of the true 
satellites. For the purposes of this test ( N int ) is estimated 
using the random points method (i.e., assuming uniform 
projected density of interlopers). 

The function w(R) can be compared to the surface 
density of the true satellite galaxies, Etruesat(-R), esti¬ 
mated using the random points method. If the estimate 
of Struesat (-R) is affected by incompleteness, the func¬ 
tions w(R) and £truesat(7?) should have different shapes 
because w(R) is defined as a ratio of quantities equally 
affected by area incompleteness, which should cancel out 
the effect. Etruesat(-R), on the other hand, will be af¬ 
fected. Conversely, w(R) and Etruesat(-R) should have the 
same shape if effects of incompleteness on Etruesat {R) are 
negligible. This is because for the random points inter¬ 
loper Subtraction: Etruesat = ( {Ns at) / {n t) l)£int = 
u;Ei n t, where E; nt (i?) = const. 

Figure 0 shows the ratio of the surface density pro¬ 
file of satellites estimated for our volume-limited sam¬ 
ple using the random points method and function w(R ), 
computed for the same sample. Over the projected radii 
test, the ratio is consistent with a constant. In addi¬ 
tion, the marginalized best-fit slopes of E trU esat and w(R) 
agree within statistical errors: a = —1.33 ± 0.09 and 
a = —1.41 ± 0.06, respectively. We therefore conclude 
that our measurements of the surface density profiles 
of satellites are not significantly affected by area incom¬ 
pleteness. 


4.3. Flux-limited samples 


TABLE 6 

Estimated Power-Law Slopes for the Flux-Limited Sample 


Input Data 

Sv = 500 km s -i 




satellite sample 


-1.27 

±0.06 

clustered random points 


-1.52 

±0.07 

nearby points (Aj? C orr = 1 - 

2) 

-1.64 

±0.07 

Sv = 1000 km s —1 




satellite sample 


-1.17 

±0.06 

clustered random points 


-1.56 

±0.07 

nearby points ( AR corr = 1 - 

2) 

-1.66 

±0.08 



Satellites per primary 



M r 
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Fig. 9. — Statistics of the satellite samples for a flux-limited 
sample. The solid line shows the results for the Sv = 500 km s -1 
criterion, while the dotted line shows the Sv = 1000 km s~ 1 cri¬ 
terion. Top: The number of satellites per primary for primaries 
with at least one satellite. Center: The r-band magnitude his¬ 
togram for satellites. Bottom: The magnitude differences between 
satellites and primaries. 


We could significantly increase the number of objects 
in our samples by eliminating the limiting absolute mag¬ 
nitude and increasing the maximum depth of the sample. 
Increasing the depth, however, will increase the minimum 
separation of the sample as dictated by the fiber collision 
problem. We, then, only eliminate the magnitude limit 
to create a flux-limited sample. For these flux-limited 


samples, we apply the same isolation and satellite crite¬ 
ria as the volume-limited sample. For the satellite sam¬ 
ple with Sv = 500 km s -1 , we have 678 galaxies. For the 
larger Sv = 1000 km s _1 sample, with a larger minimum 
separation and limiting magnitude, there are 786 galax¬ 
ies. The flux-limited samples are summarized in Table 
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R (h - 1 Mpc) 


intervals. 

The marginalized best-fit power-law slope values in Ta¬ 
ble El show similar results to those found in the volume- 
limited sample. The satellite sample with Sv = 500 km 
s _1 has a somewhat steeper slope than the Sv = 1000 
km s -1 sample, —1.27 ± 0.06 to —1.17 ± 0.06. In ad¬ 
dition, the clustered random points and nearby points 
methods produce very consistent results in both sam¬ 
ples, a = —1.52 ± 0.07 and a = —1.64 ± 0.07 for the Sv 
= 500 km s -1 sample. 

The similarity of the results for the volume-limited and 
flux-limited samples suggest that the flux-limited data do 
not induce any bias to the sample while increasing the 
statistical significance of the results. 

4.4. Trends with luminosity 



Fig. 10.— The satellite sample compared to interloper subtracted 
samples for the flux-limited sample with velocity criterion Sv = 
500 km s -1 . Top: The best-fit power-law for the satellite sample 
(solid) and the interloper subtracted profile for the ARcott = 1 - 
2 h ~ 1 Mpc nearby points (dotted). The data is shown with error 
bars: satellite sample (square points, thick Poisson error bars) and 
nearby points (open circles, thin error bars). Bottom: The 68%, 
90%, and 99% confidence intervals for the satellite sample (solid 
lines) and the nearby points sample (dotted lines). The contours 
are labeled by the unmarginalized best-fit slope of each profile. 
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Fig. 11.— The left-hand panels show the results of the nearby 
points method for samples split by bright and faint primaries. The 
right-hand panels shows the results when samples are split into 
bright and faint satellites. In all cases, the ’bright’ data is shown 
in square points, with thick error bars and the ’faint’ data is shown 
in open circles with thin error bars. In addition, the unmarginalized 
best-fit power-law for the ’bright’ data is shown with a solid line 
and a dotted line for the ’faint’ data. The top panels show results 
from the volume-limited sample and the bottom panels show the 
results from the flux-limited sample. 


01 Statistics of the satellite samples in the flux-limited 
set are plotted in Fig. 01 while the satellite sample and 
the nearby points interloper subtracted results with Sv 
= 500 km s -1 are plotted in Figure ITU1 for bins of 35 h _1 
kpc starting at the minimum separation of 32.9 h~ l kpc. 
In addition, the marginalized best-fit slopes are listed in 
Table El 

The flux-limited results are consistent with those of 
the volume-limited sample. The clustered random points 
and nearby points methods for interloper removal once 
again show steeper profiles than the satellite sample. 
We also find best-fit power-laws that are similar to the 
volume-limited sample for the satellite sample (a = 
— 1.27) and nearby points method (a = —1.64) in the 
top panel of Fig. EH The bottom panel shows the confi¬ 
dence regions for the two fits, where the slopes of the two 
distributions do not overlap within the 99% confidence 


We now test the projected radial distribution for dif¬ 
ferences between bright and faint primary galaxies and 
bright and faint satellites (see Figure ITO for the nearby 
points method and Sv = 500 km s -1 . We first split 
the primary sample into two, a bright sample of —23 < 
M r < —21 with 125 primaries and a faint sample of 
—21 < M r < —20 with 746 primaries. The volume- 
limited set has 161 possible satellites in the bright pri¬ 
maries sample and 171 in the faint primaries sample. The 
flux-limited set has 225 possible satellites in the bright 
primaries sample and 297 in the faint primaries sample. 
The left-hand panels of Figure HD show that, in all cases, 
the amplitude of the projected number density profile is 
always larger for the bright primary samples. The bright 
and faint primary samples show consistent slopes. The 
slopes marginalized over the normalization for the bright 
volume- and flux-limited samples are —1.56 ± 0.14 and 
—1.55±0.11, respectively, while the corresponding slopes 
for the faint samples are —1.57 ± 0.17 and —1.68 ± 0.10. 
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We thus do not find a significant dependence of the con¬ 
centration of the radial profiles on the luminosity of the 
primary galaxies in the luminosity range probed. 

Next, we split the satellite sample into samples of 
bright and faint satellites at the median absolute mag¬ 
nitude of each satellite sample. The volume-limited set 
is split at M r = -18.28, while the flux-limited at M r = 
-17.76. The luminosity range of satellite galaxies probed 
is not large, with the medians of the satellite samples 
~ 0.1 T*. The right-hand panels of Fig. 1 11 1 show that the 
best-fit power-law slopes, marginalized over the normal¬ 
ization, for the bright and faint satellites are consistent: 
— 1.72 ± 0.16 and —1.45 ± 0.16, respectively, for the the 
volume-limited set, and —1.55 ± 0.11 and 1.72 ± 0.11 for 
the flux-limited set. 

4.5. Putting everything together: comparisons with 
ACDM expectations 

In this section we compare the observed projected den¬ 
sity profile of satellites in our SDSS sample to mock satel¬ 
lite samples derived from high-resolution simulations of 
the concordance ACDM model. For this comparison, we 
compare to catalogs that are representative of subhalo, 
satellite galaxy, and dark matter distributions, interest¬ 
ing reference points to compare to the observational re¬ 
sults. In the first catalog, we use samples of isolated 
halos and subhalos selected simply using their maximum 
circular velocity, V max . The only criterion we use for 
subhalos is I4iax > 100 kms -1 to avoid possible biases 
due to effects of resolution for smaller subhalos. We use 
two additional samples to probe the subhalo distribu¬ 
tion given two different ways of selecting subhalos. In 
the first, galaxies and satellite galaxies are associated 
with halos and subhalos, and selection is based on the 
r-band galaxy luminosities assigned in such a way that 
the observed luminosity function of galaxies, their clus¬ 
tering, and gal axy-mass correlations a re well reproduced 
(see § |2 and iTasitsiomi et alJ 12004 for details). The 
selection criteria of the primary galaxies and satellites 
mimic those applied to the SDSS samples. This sample 
is analogous to Test Sam ple 3 ( see §13, except that_our 
sample is taken at z = 0. IConrov. Wechsler. fc Kravtsovl 
: 2005b : present a very promising modification to the lu¬ 
minosity assignment scheme discussed which, coupled 
with high-resolution dissipationless simulations, repro¬ 
duces the small-scale galaxy clustering and its luminos¬ 
ity dependence observed in SDSS remarkably well. This 
scheme assigns r-band galaxy luminosities for subhalos 
not on the D m ax at 2 = 0 but using the Vj nax at the 
epoch of accretion for the subhalo. Given that the bary- 
onic components of subhalos should be more resistant 
to the physical processes that evolve the properties of 
subhalos as they fall into their parent halo, this sample 
should better reflect the distribution of satellite galaxies. 
Finally, we use the same primaries as for the subhalo 
sample, but the satellite sample is now constructed by 
randomly sampling DM particles surrounding primaries. 

Figure lT^l shows the interloper-corrected profile for our 
volume-limited SDSS sample along with the profiles for 
the true satellites distribution for satellites assigned lumi¬ 
nosities by I4iax at £ = 0 and V max at accretion and the 
sample of dark matter particles around halos with Vmax 
values in the range 200 — 300 kms -1 . The normalization 
of the dark matter particle profile is set arbitrarily in the 



Fig. 12. — Comparison of the projected number density of satel¬ 
lites in the volume-limited SDSS sample ( open circles with error 
bars connected by the solid line) and in the three mock satellite 
simulation-derived samples (see text for details). The dot-dashed 
line represents dark matter distribution around primary halos in 
simulations, while dotted and dashed lines show satellite profiles 
using two different ways of assigning luminosities to subhalos (see 
text). 


plot. The figure shows that the shape of the profiles for 
all three simulation-derived samples is similar to that of 
the SDSS sample. The distribution for satellites assigned 
luminosities by Vmax at z = 0, however, is shallower than 
the observed distribution. The distribution for satellites 
assigned luminosities by Vmax at accretion also appears 
somewhat shallower, but less so than for the other lu¬ 
minosity assignment. The dark matter particles, on the 
other hand, have a somewhat steeper radial profile than 
the observed satellites. 

FigureEHconrpares the best fit slopes for the full range 
of radii for the different simulation-derived samples dis¬ 
cussed above with the corresponding best fit slope for 
the radial profile of the SDSS satellites. We present ra¬ 
dial distributions of subhalos and dark matter particles 
for the host halos in three ranges of V max to illustrate 
the dependence of the satellite distribution of the par¬ 
ent system mass. For the Vmax-selected subhalo catalog 
for primary systems with Vmax = 200 — 300 kms -1 , the 
slope and normalization of the distribution are in reason¬ 
able agreement with the data at the 2cr level. The slope 
is somewhat shallower (a ss —1.3 — 1.4) than observed, 
but overlap with the latter within the uncertainties. The 
slope of the corresponding catalog in which satellites are 
selected from DM particles, exhibit a somewhat steeper 
slope (a ~ —1.8) than observed for the SDSS galaxies, al¬ 
though the difference is not statistically significant. Note 
that the normalization, A , for this sample is not mean¬ 
ingful, while it is meaningful for the sublralo sample. For 
primary systems with larger Vmax ranges, the subhalo ra¬ 
dial disribution becomes progressively shallower than the 
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Fig. 13. — The 95.4% confidence intervals for the best-fit power- 
laws of the interloper subtracted data and for the true satellite dis¬ 
tributions of several samples drawn from the simulations. The two 
filled contours represent the volume-limited sample (on the left) 
and the flux-limited sample (on the right). The solid and dashed 
contours show satellite samples using two different luminosity as¬ 
signment schemes, as discussed in the text. Samples are drawn 
from Vmax selected primary samples for subhalos with V m &x >100 
km s -1 (dotted contour lines) and for DM particles (crosses), with 
the Vmax ranges of 200-300, 300-400, and 400-600 km s —1 showing 
increasingly shallower slopes and larger amplitudes. The normal¬ 
ization, A, for samples of DM particles is not meaningful, as the 
samples are drawn from arbitrary fractions of DM particles around 
host halos. 


data because the fixed projected radial range used probes 
different radii with respect to the halo virial radii. Most 
of the primaries in our sample should correspond to ha¬ 
los with 14iax < 300 bis -1 , so it is not surprising that 
the data matches the subhalo distribution for primaries 
with Vmax = 200 — 300 kms -1 the best. For primaries 
of this size, the virial radius is expected to significantly 
smaller than 500 h -1 kpc, the maximum projected radius 
measured. Simulations, however, show bound subhalos 
at several times the virial radii of the host halos. 

Figure EH also compares the best fit slopes for the 
SDSS to satellite samples with different luminosity as¬ 
signments. Both assignments are in agreement with the 
best fit slope of the SDSS volume-limited sample at the 
2cr level. Assigning luminosities based on a subhalo mass 
at accretion, however, results in a steeper slope than as¬ 
signing luminosities based on a subhalo mass at z = 0 
and is more consistent with the data. 

As previously discussed, the slopes of satellite surface 
density profiles around SDSS galaxies are in between 
those measured for subhalos and those for dark matter 
distributions in simulations of the concordance cosmol¬ 
ogy, and the measured slopes are closer to those mea¬ 
sured for dark matter distribution. In addition, the bias 
in the interloper subtraction methods suggest that the 
true satellite distribution is somewhat steeper than our 
best fit value and even more compatible with the dark 
matter distribution. Our results are in qualitative agree¬ 
ment with the recent observational studies on the ra¬ 
dial distribution of galaxies in groups and clusters (e.g., 


TABLE 7 

Estimated Power-Law Slopes for Different Samples 


of True Satellites 

Primary+satellite sample R m in < R < 0.5 

0.1 < R < 0.5 

SDSS (volume) 

SDSS (flux) 

-1.58 ±0.11 
-1.64 ±0.07 

-1.77 ±0.23 
-1.72 ±0.17 

Mr (Vmax a.t Z = 0) 

M r (Vinax Ut Z a cc ) 

Vmax = 200-300 (DM) 

Vmax = 200-300 (subs) 

Vmax = 300-400 (DM) 

Vmax = 300-400 (subs) 

Vmax = 400-600 (DM) 

Vmax = 400-600 (subs) 

-1.12 ±0.10 
-1.31 ±0.09 
1.895 ±0.010 
-1.40 ±0.07 
1.705 ±0.009 
-1.23 ±0.06 
1.513 ±0.008 
-0.94 ±0.04 

-1.36 ±0.18 
-1.38 ±0.17 
-2.111 ±0.020 
-1.65 ±0.13 
-1.864 ±0.018 
-1.65 ±0.11 
-1.674 ±0.015 
-1.26 ±0.07 


Note - Samples using isolated halos chosen by Vmax and sub¬ 
halos assigned a r-band luminosity use binning of 70/i _1 kpc, 
with first bin starting at _R m in = 25.6/i 1 kpc. Subha¬ 
los in these samples are selected to have circular velocities 
Vmax > 100 km s _1 . Data samples use binning previously 
described in the text. 


Lin et al. 200J 

Hansen et al. 2005: Collister A/; Lahav 

20051 Yang et a,l. 

2005; Coil et al. 2006), which find con- 


centrations of galaxy radial profiles somewhat lower than 
the concentrations expected for the matter distribution 
of their parent halos. However, here we seem to find 
that the difference between the DM distribution and the 
satellite distribution is not as large in galaxy-sized sys¬ 
tems compared to groups and clusters. 

Several previous studies which considered satellite 
distribution in cluster-sized, halos showed that sub¬ 
halos appear to have more extended and shallower 
radial distributions c ompared to that o f dark matter 
|_|1998|: |Colfnet, afl_|1999|_jGhigna et^lj 

2000[__|S]ji i iiiKde^a] |__|2001l iDeLuci^et^LLJ^OO^ 

DigmantLst al .112004 i Gao et alJl2004lNagai fe Kravtsov! 

20051) . Tidal evolution and merging modify the subhalo 
profile, especially within inner ~ 50% of the virial 
radius, primarily because they modify properties of 
subhalos, such as its bo und mass or circular velocity 
llNagai fc Kravtsov! 12005 1. For galaxy-sized systems in 
our simulations, then, we would thus expect subhalo 
profile to be significantly flattened at R < 100h -1 kpc. 
Table □ presents the best-fit slopes for the full range 
of radii and for the fits restricted to radii between 0.1 
and 0.5 h~ x Mpc for all the data and simulation samples 
discussed. The table shows that the best-fit power-laws 
for bins of R > 100/i _1 kpc are steeper than for those 
including smaller projected radii. This is generically 
true, but in particular subhalo distributions seem to 
steepen more than the DM distributions, suggesting 
that subhalo distributions are flattened within ~50% of 
the virial radius as in cluster-sized halos. For example, 
for subhalos assigned luminosities at z = 0, the slope for 
the restricted set of radii is —1.36 ± 0.18 - « 0.2 steeper 
than the fit using all data points - and the slope steepens 
further if we constrain the fit to even larger radii, where 
the best fit slope values of DM and subhalo distributions 
agree within error bars (note, however, that at these 
radii the fit errors become larger). On the other hand, 
the slope for subhalos assigned luminosities at the epoch 
of accretion and the observational results steepen less, as 
would be expected given that the baryonic components 
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of subhalos should be more resistant tidal evolution 
and merging. We attribute some of the differences in 
slopes and shallower slopes for our simulation satellite 
samples to tidal evolution effects, and to the differences 
in object selection in s imulation and SDSS samples 
llNagai & Kravtsov! 12005 1. A better understanding of 
the differences between selection criteria is necessary for 
more stringent tests of the theoretical predictions. 

5. DISCUSSION AND CONCLUSIONS 

Modern large galaxy redshift surveys allow one to 
study the distribution of satellites around galaxies and 
clusters with unprecedented statistical power, while con¬ 
trolling biases and completeness in a systematic way. 
In addition, redshift information can be used to select 
galaxies from relatively non-crowded environments and 
to account for interloper contamination in a rigorous way. 
Cosmological simulations are also sufficiently mature and 
allow systematic tests of the interloper subtraction algo¬ 
rithms. Galaxy and satellite samples, for instance, can 
be constructed to mimic observational selection criteria. 
Examples of studies using such surveys are constraints 
on ^the JDM^ halos ^of galaxies ^fro m satellite kinematics 

and the anisotropy of the distribution of satellite galaxies 
llSales fc LambasH2004 lBrainerdll2005al) . 

In this work, we use the SDSS spectroscopic survey 
to estimate the projected radial distribution of satellites 
around isolated primaries. We use areas of the survey 
which are at least 90% complete and check for the ef¬ 
fects of incompleteness by comparing the surface density 
profile and conclude that our results are not affected by 
the (small) incompleteness of the sample. We construct 
samples of primary and satellit e galaxies with isol ation 
criteria similar to those used bv lPrada et al . (l2003t) . We 
use high-resolution cosmological simulations of the con¬ 
cordance ACDM cosmology to develop and carefully test 
new methods of correcting for interloper contamination. 
Our main results and conclusions can be summarized as 
follows. 

1. Using mock galaxy catalogs derived from high- 
resolution cosmological simulations, we show that 
interlopers can significantly bias the shape of the 
projected surface density profile of faint satellites 
around bright galaxies, making it shallower (bi¬ 
asing the power-law slope a of the radial profile, 
T,(R) oc I?“, by Aabias > 0.5). We also show that 
the most straightforward methods do not correct 
interloper contamination properly. For example, 
the random points method, which assumes uniform 
distribution of interlopers in space, underestimates 
the fraction of interlopers in the satellite sample 
by oversampling voids compared to clustered areas 
where most galaxies in the sample reside. 

2. We develop two new methods to account for the 
interloper contamination: the clustered random 
points method and the nearby points method, vari¬ 
ants of the random points method, designed to 
sample environments similar to those of the clus¬ 
tered galaxies in the observed samples. Tests on 
the mock samples show that the methods perform 
consistently well, reducing the interloper bias in the 


best-fit power-law slopes of the satellite profiles to 
only Aabias ~ 0.1 for the nearby points method. 

3. We apply these methods in our analyses of the 
volume- and flux-limited SDSS spectroscopic sam¬ 
ples. The best fit power-law slope for the volume- 
limited SDSS satellite sample, after interloper con¬ 
tamination correction, is —1.58 ±0.11 in the range 
of projected separations of 32.9 < R < 500 h~ l kpc. 
Note that we estimated a systematic bias in the 
derived slope of Aabias ~ 0.1 for nearby points 
interloper subtraction method, which implies that 
the true slope of the SDSS satellites may be a = 
a raw — Aabias ~ —1.7. We find similar values of the 
best fit slopes for the flux-limited samples, and for 
samples of primary galaxies in different absolute 
magnitude ranges. We thus do not find evidence 
for the dependence of the shape of satellite radial 
distribution on the luminosity of their host galaxy. 

4. Comparison of the observed radial distribution of 
the SDSS satellites to the distribution of sublralos 
and dark matter around galactic halos in dissipa¬ 
tionless ACDM simulations shows that the slope 
of the SDSS satellite radial profile is in between 
those measured for subhalos and for dark matter 
(closer to dark matter). Subhalos thus appear to 
have more extended and shallower radial distribu¬ 
tions than the observed satellites. The dark matter 
distribution is somewhat steeper than the observed 
satellite profiles, but the difference is not statisti¬ 
cally significant. 

Recently, Ivan den Bosch et _yj. 11209511 and 
iSales fe Lambadi (I2005H studied the projected radial 
distribution of satellite galaxies around isolated galaxies 
using the Two D egree Field Ga la xy R edshift Survey 
(2dFGRS). Ivan den Bosch et alJ ( I2005D used mock 
galaxy redshift samples derived from large cosmological 
simulations to develop an iterative method of interloper 
rejection for the 2dFGRS and found that the data is 
generally consistent with the dark matter profile at 
large projected radii, but concluded that incompleteness 
of close pairs in the su rvey prevent strong constraints. 
ISales & Lambadl (12005}) used isolation criteria for the 
primary galaxies, enforcing that any neighbor within a 
region of 700 h~ l kpc and AU = 1000 km s -1 must be 
at least 1 magnitude fainter. Satellites were assumed to 
be any object within projected separation of 500 kpc and 
AU = 500 kms^ 1 that is at least 2 magnitu des fainter 
than the host galaxy. ISales fc Lambas ; (l2f)05t) estimated 
the effect of close-pair incompleteness by deriving a 
control sample of projected satellites with velocity 
difference of 2000 < AU < 10,000 kms" 1 . At R > 20 
hr 1 kpc, where the number density profile of the control 
sample is flat, the satellite s amples were con s idered 
complete. For the total sample, ISales fc Lambas ( 12005} ) 
quoted the best-fit power-law of slope a = —0.96 ± 0.03. 
This value, however, was derived without applying 
correction for interloper contamination. Note that we 
derive signi ficantly steeper p r ofiles of satellites compared 
to those of ISales fc Lambadi feOOSl) . which we attribute 
to the rigorous correction for interlopers and to possible 
fiber-collision bias in the 2dFGRS. 
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The results presented in this study provide interest¬ 
ing hints of the possible differences between observed 
satellite distributions and the expected distribution of 
subhalos in their parent halos. However, the statistical 
errors are still rather large. Significant improvements in 
the statistics are needed to address this question further. 
Larger statistics would also allow us to go beyond the 
average profiles and study the distribution of satellites 
as a function of satellite (e.g., luminosity and color) and 
host galaxy properties. 
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